Identified Particle Spectra and Anisotropic Flow in an Event-by-Event Hybrid 
Approach in Pb+Pb collisions at ^/snn = 2.76 TeV 
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The first results from heavy ion collisions at the Large Hadron Collider for charged particle spectra 
and elliptic flow are compared to an event-by-event hybrid approach with an ideal hydrodynamic 
expansion. This approach has been shown to successfully describe bulk observables at RHIC. With- 
out changing any parameters of the calculation the same approach is applied to Pb+Pb collisions 
at ^/snn = 2.76 TeV. This is an important test if the established understanding of the dynamics of 
relativistic heavy ion collisions is also applicable at even higher energies. Specifically, we employ the 
hybrid approach with two different equations of state and the pure hadronic transport approach to 
indicate sensitivities to finite viscosity. The centrality dependence of the charged hadron multiplic- 
ity, pt spectra and differential elliptic flow are shown to be in reasonable agreement with the ALICE 
data. Furthermore, we make predictions for the transverse mass spectra of identified particles and 
triangular flow. The eccentricities and their fluctuations are found to be surprisingly similar to the 
ones at lower energies and therefore also the triangular flow results are very similar. Any devia- 
tions from these predictions will indicate the need for new physics mechanisms responsible for the 
dynamics of heavy ion collisions. 

PACS numbers: 25.75.-q,24.10.Lx,24.10.Nz 
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Recently, the first results from heavy ion collisions at 
E cm = 2.7 6 A TeV have been published by the ALICE 
collaboration [3-[l] ■ To study strongly interacting matter 
at high temperatures has been the goal of the relativistic 
heavy ion program at the Relativistic Heavy Ion Collider 
(RHIC) since more than a decade. The 10 times higher 
beam energies at the Large Hadron Collider (LHC) al- 
low for the investigation of the dynamical evolution of 
nucleus-nucleus collisions that have been established in 
Au+Au collisions at E cm = 200 A GeV in a different kine- 
matic range Q . It is especially interesting, if the matter 
still behaves as a almost perfect liquid or if the quark 
gluon plasma becomes more viscous by going to higher 
temperatures [6|-|9j. 

Hybrid approaches that are based on hydrodynam- 
ics for the hot and dense stage of the evolution that 
is coupled to hadronic transport approaches to describe 
the successive decoupling of the matter have been very 
successful in describing the properties of the bulk mat- 
ter that is created at RHIC (TMl and LHC [H|. 
Within the last year, full event- by- event hydrodynamic 
approaches [Pa - l21j have become more favorable because 
it has turned out that the effect of initial state fluctua- 
tions on final flow observables needs to be studied in a 
consistent way to draw quantitative conclusions for e.g 
the value of the shear viscosity. 

It is important to apply well-established approaches for 
the dynamical evolution of heavy ion reactions at RHIC 
to the higher energy collisions at LHC without tuning pa- 
rameters to investigate how good the energy dependence 
is described by a specific model. Differences between the 
experimental data and the theoretical calculations im- 
ply the need for new physics concepts to be applied. To 
understand the bulk evolution at LHC energies is a pre- 
requisite for all further detailed studies of e.g. jetquench- 
ing and energy loss or electromagnetic probes [22[ HH ■ 



In this manuscript a event- by- event hybrid approach 
based on the Ultra-relativistic Quantum Molecular Dy- 
namics (UrQMD) transport approach with an embedded 
(3+1) dimensional ideal hydrodynamic evolution is ap- 
plied to lead-lead collisions at LHC energies. First, a 
comparison to the available experimental data on charged 
particle multiplicities, px spectra and elliptic flow is car- 
ried out with the exact same parameter set that has been 
applied at RHIC energies. Then, transverse mass spectra 
and triangular flow for identified particles are predicted. 
To compare the amount of initial state fluctuations, the 
probability distributions of coordinate space eccentrici- 
ties are compared to their corresponding values at RHIC. 

Let us start with a short description of the hybrid ap- 
proach that has been developed for SPS energies [l8[ and 
recently successfully applied to gold-gold collisions at the 
highest RHIC energy [2J-|26[ . The early non-equilibrium 
evolution is described by the UrQMD approach [27], HH , 
where the two lead nuclei are initialized according to 
Wood-Saxon profiles followed by binary interactions of 
the nucleons. The main contribution to the particle pro- 
duction at high energies is achieved by string excitation 
and fragmentation where for the hard collisions (momen- 
tum transfer Q > 1.5 GeV) are treated by PYTHIA [H- 

& 

At the so called starting time of i s tart = 0.5 fm the 
particle distributions are transferred to energy, momen- 
tum and net baryon density distributions by representing 
each particle with a three-dimensional Gaussian distribu- 
tion (width (7 = 1 fm) that is Lorentz-contracted along 
the beam direction. During the following ideal relativis- 
tic one fluid evolution [32|, [33| two different equations of 
state are employed, one representing a hadron gas (HG- 
EoS) and one including a cross-over deconfinement phase 
transition based on a chiral approach (DE-EoS) [3J, [35[ . 
The transition back to the hadronic transport approach 
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happens on a constant proper time hypersurface, where 
the Cooper-Frye equation is applied on transverse slices 
of thickness Az = 0.1 — 0.2 fm that have cooled down be- 
low an energy density of 5erj ~ 730 MeV/fm 3 [36]. This 
approach provides the full final phase space distributions 
of the produced particles for each event and can be com- 
pared to the pure transport approach by turning off the 
hydrodynamic evolution which allows for a qualitative 
study of viscous effects. 
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FIG. 2: (Color online) Transverse momentum spectra of 
charged particles for four different centralities calculated in 
the UrQMD transport and the hybrid approach compared to 
the available experimental data [2|. 



FIG. 1: (Color online) Charged particle multiplicity at midra- 
pidity (\rj\ < 0.5) as a function of the number of participants 
in Pb+Pb collisions at ^/snn = 2.76 TeV calculated in the 
UrQMD transport and the hybrid approach compared to the 
experimental data pj|. 



Centrality class 


Impact parameter range 


0-5% 
5-10% 
10-20% 
20-40% 


b < 3 fm 
b = 3 - 5 f m 
b = 5 - 7 f m 
b = 7 - 10 fm 



The first observable to look at is the charged particle 
multiplicity at midrapidity. In Fig. [TJ the calculation of 
the centrality dependent multiplicity scaled by the num- 
ber of participants (estimated in a Glauber approach) is 
shown. The hadronic transport approach UrQMD pro- 
vides a reasonable description of the multiplicity. For 
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central collisions the predictions published in 
right on top of the ALICE data while with decreasing 
centrality the number of charged particles is a little lower 
than in the data. This fair agreement with the data hints 
to the fact that the main particle production can be de- 
scribed by the initial binary nucleon-nucleon interactions 
treated by PYTHIA. The hydrodynamic evolution does 
not affect the particle production. Since ideal hydrody- 
namics implies an isentropic expansion this means that 
the charged particle multiplicity is determined in the ini- 
tial state and by the final resonance decays. 

For the following calculations of spectra and collective 
flow four different centrality classes have been chosen that 
match the ones applied by the ALICE collaboration as 
they are listed in the following table: 



The transverse momentum spectrum for charged par- 
ticles in the mentioned centrality classes are compared to 
experimental data in the most central bin (see Fig. 
The main difference between the hybrid and the trans- 
port calculation is in the slopes of the spectra. As ex- 
pected the hydrodynamic evolution leads to a purely ex- 
ponential pt dependence which describes the data until 
Pt < 3 GeV very well. At higher transverse momenta 
the power law tail from hard processes becomes impor- 
tant for a good agreement with the measured values. In 
the range from 4 to 6 GeV the non-equilibrium descrip- 
tion exemplified by the UrQMD calculation provides a 
better description of the experimental data. 

In Fig. |3] predictions for the transverse mass spec- 
tra at midrapidity of pions, kaons and protons are pre- 
sented. The pion spectra are very similar to the charged 
particle spectra since they represent the major fraction 
of the newly produced particles in the collision. Kaons 
are strange mesons and protons are chosen because they 
have a higher mass and are baryonic degrees of freedom. 
The general features of the transverse mass spectra are 
similar to the ones observed at RHIC and imply a col- 
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FIG. 3: (Color online) Transverse mass spectra of negative 
pions (top), positive kaons (middle) and protons (bottom) for 
four different centralities calculated in the hybrid approach 
with two different equations of state. 



lective radial velocity that drives all the particle species. 
The two different equations of state lead to very similar 
results with the deconfinement transition having a little 
steeper slope due to the more rapid expansion due to the 
higher pressure in the quark gluon plasma phase. 

After proving a rather successful agreement with basic 
quantities like the multiplicity and transverse momen- 
tum spectrum the next step is to look at anisotropic 
flow observables. The elliptic flow has been calculated 
with respect to the reaction plane by averaging over all 
charged particles in all events to be compared to the 
ALICE measurement that relies on the four-particle cu- 
mulant method in two centrality bins. Fig. [4] shows a 
good agreement between the hybrid calculations and the 
data, especially between pt=0.8-2.5 GeV. In the very 
low transverse momentum region the hybrid approach 
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FIG. 4: (Color online) Elliptic flow of charged particles as a 
function of transverse momentum for four different centralities 
calculated in the hybrid approach with two different equations 
of state compared to the experimental data[||. 



underpredicts the data which has been observed in other 
calculations as well [l4| • At higher px again the influence 
of hard processes needs to be taken into account. 

To quantify the shape of the initial conditions em- 
ployed for the hydrodynamic calculation and its event- 
by-event fluctuations Fig. [5] shows the probability distri- 
bution of the coordinate space asymmetry characterized 
by the eccentricity and the triangularity as defined in 
j25| . The initial e n coefficients have been calculated in 
each event and the normalized probability distribution is 
plotted for two different centrality bins. 

For central collisions the mean value and the shape of 
the distributions are very similar for the participant ec- 
centricity and the triangularity since both of them are 
mainly generated by fluctuations. For more peripheral 
collisions the eccentricity is influenced by a large geome- 
try component due to the ellipsoidal shape of the initial 
state in the transverse plane. Therefore, the mean ec- 
centricity is larger and the fluctuations increase leading 
to a wider distribution, whereas the triangularity stays 
smaller and the distribution has a smaller width. 

Since the triangularity has been introduced because of 
its sensitivity to initial state fluctuations the higher mul- 
tiplicity at LHC energies triggers the expectations that 
the fluctuations become smaller compared to RHIC en- 
ergies. In Fig. [5] the triangles and diamonds depict the 
eccentricity and triangularity calculation from UrQMD 
initial conditions for Au+Au collisions at E om = 200A 
GeV. Surprisingly, the e„ distributions match almost ex- 
actly the ones at LHC energies for the two similar cen- 
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FIG. 5: (Color online) Probability distributions of the ec- 
centricity coefficients (62/63) in the UrQMD initial state. The 
diamonds and triangles indicate the corresponding distribu- 
tions for Au+Au collisions at RHIC at v / %n = 200 GeV. 



perimental result for minimum bias Pb+Pb collisions at 
E cm — 2.76A TeV. Both elliptic and triangular flow show 
the expected mass splitting between pions and protons. 
Furthermore, the elliptic flow has roughly double the 
value of the triangular flow result as it has previously 
been calculated for RHIC energies. The anisotropic flow 
are not very sensitive to the choice of the equation of 
state which points to the fact that the interplay of freeze- 
out transition and flow development due to pressure gra- 
dients results in very similar values independent of the 
existence of a quark gluon plasma phase. 
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FIG. 6: (Color online) Triangular flow of pions and protons 
as a function of transverse momentum compared to the corre- 
sponding elliptic flow result calculated in the hybrid approach 
with two different equations of state for minimum bias Pb+Pb 
collisions at Jsnn = 2.76 TeV. 



To calculate the transverse momentum dependence of 
the anisotropic flow coefficients for identified particles 
(see Fig. [6]) the event plane method including removal 
of auto-correlations and resolution correction has been 
applied as described in more detail in [25j . In this way, 
the predictions can be directly compared to a future ex- 



FIG. 7: (Color online) Triangular flow of charged particles as 
a function of transverse momentum calculated in the hybrid 
approach for four different centrality classes of Pb+Pb colli- 
sions at ^/snn = 2.76 TeV compared in two centrality bins to 
the corresponding results of the UrQMD transport approach. 

Finally in Fig. [7]predictions for the transverse momen- 
tum dependence of triangular flow for charged particles 
in the four different centrality classes are made. There 
is a weak centrality dependence, but overall the differ- 
ent curves from the hybrid calculation look very similar 
qualitatively. Since the calculation with the deconfinc- 
ment phase transition leads to very similar results, only 
the hadron gas result is shown in this figure to reduce 
the number of lines displayed. To investigate the effect 
of a finite viscosity during the expansion the calculation 
is compared to the pure UrQMD result. Within the very 
viscous hadronic transport approach the pressure gradi- 
ents are too small to transfer the initial coordinate space 
asymmetry to a final state momentum space asymme- 
try. For impact parameters up to b = 5 fm the result is 
consistent with zero and even for very peripheral event 
the result is negligible especially in the low transverse 
momentum region where most of the particles are. This 
leads to the conclusion that a finite triangular flow mea- 
surement is a strong indication for almost ideal hydro- 
dynamic behavior during the hot and dense evolution of 
the heavy ion reaction. 

To summarize we have presented a comparison of bulk 
observables measured in Pb+Pb collisions at E cm = 
2.76A TeV to a state-of-the-art event- by- event hybrid 
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description. Employing the same parameters as have 
been used for Au+Au collisions at RHIC a reasonably 
good agreement for the multiplicity, transverse momen- 
tum spectra and elliptic flow is achieved. One can con- 
clude that the basic assumptions about the major ingre- 
dients that are needed to describe the dynamic evolution 
of heavy ion reactions do not have to be revised at LHC 
energies. 

With a non-equilibrium initial state, an ideal hydrody- 
namic evolution and a hadronic afterburner predictions 
for identified particle transverse mass spectra and elliptic 
as well as triangular flow are made. The initial state fluc- 
tuations quantified by the probability distribution of the 
initial eccentricity and triangularity are very similar at 
LHC and RHIC energies within the UrQMD approach. 
A large viscosity during the evolution results in almost 
negligible higher harmonic coefficients, so the triangular 
flow measurement at LHC will provide a robust confir- 



mation of the almost ideal hydrodynamic expansion. 
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